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1. INTRODUCTION 

The Unified Power Flow Controller (UPFC), is among the FACTS devices that have attracted the 
attention of many researchers because it is capable of simultaneously and independently controlling the flow 
of active and reactive powers in a network. The UPFC combines a shunt compensation, (Static Compensator 
or STATCOM), and a series compensation (Static Synchronous Series Compensator or SSSC) and has the 
ability to control three parameters associated with the transit of powers namely the line voltage, 
the impedance of the line and the load angle. The UPFC is placed on the transmission line between the source 
and the load as shown in Figure 1 [1, 2]. 





Figure 1. Electrical network with UPFC 
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It consists of two voltage source converters (VSC) connected through a common DC link. The first 
converter (A) is connected in series and the second (B) in parallel to the line (Figure 2). Converter A 
performs the main task of the UPFC by injecting an AC voltage adjustable via transformer T2. The role of 
converter B, is to supply or absorb the active power required by Converter A to the common DC circuit. 
It can also generate or absorb reactive power. 
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Figure 2. UPFC components and structure 


This paper presents the control design approach based on control Hœ and a robust controller is 
formulated in the state space domain. Hoo robust control theory is probably the control theme which has led 
to the largest number of publications and greater efforts since the mid-80 s. The main reason for its 
popularity is because it represents a very natural way to formulate the problem of robustness. 

It all started with an article by Zames [3], published in 1981 and followed by other articles [4, 5]. 
These articles did not concern the robustness problem, but rather the disturbance rejection. It is Kimura [4] 
who, in 1984, formulated the first robust control problem in terms of Hoo. Although his position of the 
problem is hardly different from that proposed by Doyle and Stein [6, 7] (novelty residing in the explicit use 
of the Hoo framework that allowed Kimura to solve the synthesis problem). Thanks to the important concept 
of the standard problem that the work of Francis Doyle and [8] were unified. 

The resolution of the standard problem grew very significantly in 1988 with Glover-Doyle 
algorithm [9], which uses state space representation. In the solutions of the problem, they introduced Ricatti 
equations [6, 8]. This approach is employed in this paper for solving the Hoo control problem. Robust 
stability is defined as the ability of a system to remain stable when subjected to perturbation such as 
modeling errors, measurement errors and external disturbances. Robust performance, on the other hand 
reflects the ability of the system to maintain the specified performance characteristics (stability, decoupling, 
time response...) when subjected to disturbances. Indeed, a physical system generally has non-linear 
characteristics that are not usually included in the model for simplification purposes. So an invariant model 
cannot accurately represent the reality and for this it is necessary to consider these errors in all control 
techniques. The Hoo control method is considered to be a very powerful design technique which has attracted 
many researchers from the electric power community. 


2. MODELING OF THE UPFC 
The simplified phase circuit of the UPFC is shown in Figure 3. Using Kirchhoffs laws, we can 


write: 
, di 

Va Va Va = Ly Hu ti) L tt idit Ri, 

dt dt C; 

dj. . di, 1... : (1) 
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L is the total inductance of the line and the load. 
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Figure 3. Equivalent circuit of the UPFC 


Using Park transformation, this system becomes: 


fi 


| \ di 
Vy Vy Vy =L, “li, +4, ae a obi, +V,, +Ri, 


E . ee di, . 
Ma Vig Va 1,7 +i toL i++ ai Peg R, 


for the shunt compensator: 
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ee : Active power consumption of the network through the parallel compensator and provided to the DC 


circuit and series compensator. 
P, : Power injected by the series compensator in the network. 


Using Park transformation: 
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Let 6 and 0 be the phase shifts between the reference and the converter output voltages V, and V,, 
respectively. The d-q components can be expressed as follows [6]: 


Vig = ky. Vde .Cos(ó); Veg = ko. Vde .sin( ô) (6) 


Via = ko. Vde . Cos( 9); Vog= Ko . Vde .sin(O) (7) 
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Assuming that the voltages at the source and receiving end are equal and the influence of the output 
shunt is neglected, the previous equations become: 














di di 
L, “+L -øL ig -OLig + eq +U cq + Rig =0 
dt dt (8) 
1 4 ^ L ij +@Li U., +Ri, =0 
—+L—+aL, ig + @Lig +Ueg t+U eg + Rig = 
S H di s'd d “cq cq q 
HO ei =0 (9) 
u 
ig -Cr —*-@C_ ucq =0 
Rearranging (8) and (9): 
dig R. 
doa er aea 
di R 
q l 
Abs oe Og ae +U cq) (10) 
d d 
Cr C AE el ig + Ocg 
duc du 1 
q f cq 
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with: Lp = L+ L,. 
This system can be written in the state space form as follows: 


E an ol, 9 -4L 
F L; Lr Lr 
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and the system transfer functions can be easily obtained from the state space equations: 


H _ pare ee 
Q G21 (s) G22 (8) || Vea 


with: 
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The instantaneous active and reactive powers generated and absorbed are defined as follows: 
3 , 3 
Fome” (Vsa lsd + Vsq sq) ; O ource = (Veq ba — Vsa ls) (1 1) 
2 2 
and the active and reactive powers absorbed by the load are: 
3 5 m 3 P X 
P,= — (Vaia + Vaia); Qr =— (V i a— Via Ug) (12) 
2 2 
with: 


' ' 


ig =iat ig and isg =ig+ Ug 


3. CONTROL OF THE UPFC 
The combination of both parallel and serial converters with a DC link provides four quadrants 
control. Figure 4 show the block diagram of overall control scheme of the UPFC [10-12]. 


Vsq 





Reference Reference 
currents currents 
generation 


generation 





7 = vV c ` 
q; O + R..||<«— a “4 Pr qr 
+ + 
Vac 


Figure 4. Full control of the UPFC 


3.1. Hoo Control Approach 
3.1.1. The Hoo Optimal Control Synthesis 

Our system is represented by a transfer matrix G(s) with a number of disturbance elements 
associated with the environment of the physical system (interference signals, etc.) and modeling errors 
(reduced order model, idealization actuators, parametric uncertainties, etc.) as shown in Figure 5 [7]. 

The aim of compensation K(s) is to ensure the stability of the closed loop system and a nominal 
satisfactory behavior. All controller qualities should be preserved as much as possible in the presence of 
external perturbations (w) and modeling uncertainties A(s), the latter being translated using interference 
signals v’. At this level, weighting matrices can be introduced on the signals (v’) and (w) to perform a 
frequency and distribution of their structural effects. Assume first that these weights have been addressed in 
the P(s) model as shown in Figure 6. 
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Figure 5. The canonical robust control problem Figure 6. A simplified representation of the control 
structure 






The goal is to find a dynamic compensator such that the stability of the system stability is ensured 
and that some transfer norm from w to z denoted F;(P, K) characterizing the performance criteria and/or 


robustness will be minimized. The problem is then: 


min | F, (P, K)| (13) 


Let f 11 2) be a partition of matrix P, the dimensions of each sub-matrix matches with the 
Pai P22 


number of inputs and outputs. Then FI(P, K) is expressed as: 
z = FI(P,K).w = (P11 + P12 KU — P22 K)-1 P21).w (14) 
The problem formulated in (35) is therefore re-written as: 
min || P11 + P12 K([— P22 K)-1 P21) || (15) 


The stabilization of the system P by the compensator K is one of the objectives of the compensation. 
However, P is a multivariable system (multiple inputs and multiple outputs) then, the concept of stability 
must be specified. 

Based on Figure 7, the standard problem is as follows: find the compensator K(s) stabilizes P(s) and 
minimizes | Tzw || 00, 

Tzw: represents the transfer matrix between w and z. 
P(s): 1s the augmented system that allows us to generate z, corresponding to the objectives of the command 
and can be represented in state space form by: 


(16) 





or by its transfer matrix. 


ple aE on ee (17) 
P, (s) P,,(s) P, (s) (Ps) 


To achieve these results, our control problem should be transformed into the standard form 
according to the principle illustrated in Figure 7. 
The objectives of the control are defined as follows: 
a) Disturbance rejection and tracking error: 
Minimize || W;(s).7(s) || œ. 
b) Noise attenuation m(s), return to maximizing margin multiplicative stability (multiplicative uncertainty 
output) 
Minimize || Ws(s).7(s) || œ. 
c) Maximize the margin of stability additive, i.e. limiting the amplitude of the control signal. 
Minimize || Wu(s).K(s). S(s) || oo. 
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This implies the minimization of: 


W, (s).S (s) 
W, (s) T(s) |œ (18) 
W (s).K(s). S(s) 


The standard problem can be represented by the equivalent diagram shown in Figure 8. 





Figure 7. Classical feedback control system Figure 8. Closed-loop system with modeling error 


which gives: 


W, —W;.G 
aa O WG ||w (19) 
y 0 W, u 
I -G 
replacing u by Ky: 
Wgs.S 

Tzw =| Wr T (20) 

W,.K.S 


The solution of this problem is obtained using the function Hj,¢ or Hing-op from the Matlab Robust 
Control Toolbox called Mixed Sensitivity Problem. Other representations of the augmented system to several 
inputs (Ww), w2, w3), and one output z can be addressed. 
The synthesis of H,, control may be summarized by: 
a) Translate the objectives oo | | i 
b) Select the weighting functions in terms of their frequency response. 
c) Set up the matrix P(s) of the equivalent standard problem. 
d) Solve the optimization problem. 
e) Test the performance of the closed-loop system. 


3.1.2. Solution of the Standard Problem and Selection of the Weighting Functions 

Doyle and others [7, 13-17] have solved this problem by performing standard programs to give the 
controller status of state space from the system and weighting functions. The assumptions used in the 
resolution of this problem are: 
1) (A, B2) stabilizable and (A, C2) detectable. 
2) Rank (D2) = m = dim (u) and rank (D21) = P2 = dim(y). 
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3) Da- | D» = [0, 7 


A-joo B, 


4) Rank 
C, D,, 


= n + m, V weļl0, of 


5) Rank|4TJeo B | -n+ P2, Yæ elo, of 
C2 D, 


a) Assumption (I) is necessary and sufficient for the existence of a solution. 

b) Assumptions (II), (IV) and (V) ensure that the problem is well posed, in other words, the 
compensator is appropriate. 

c) Assumptions (IV) and (V) indicate that Pj» (s) and P;(s) have no zeroes on the imaginary axis, 
they can be avoided [7]. 

d) Assumption (IIT) simplifies the solution. 

e) If assumption (II) is satisfied, it can used to check system (III). This operation is performed by the 
Hing program Matlab. 

f) The only problem that may arise is that of the validation of assumption (II). 


3.2. Validation of the Assumption (ID) 

If, by executing the instruction Hinf (Robust Control Toolbox of Matlab), the error signal «matrix 
Dj,» is not in full column rank», this means that the hypothesis (II) is not verified, and therefore the transfer 
functions P;2(s) and P>;(s) have zeros at infinity [9]. 

In order to ensure the closed loop stability and simultaneously achieve the desired control 
performance under process variations or in the presence of other disturbances weighting functions dependent 
of the frequency in the process are introduced as shown in Figure 9. 


Reel System 


ys) 





Figure 9. Control system with modeling error 


The controller optimization is defined by the following equation: 
|| Ws). TCs) || 0 < y Y weE[0, œ] (21) 


These weighting functions define the frequency characteristics of the signals of the system as well 
as their amplitudes. Putting |W-0@)| = 1 for a certain frequency range, the gain reduction of the 
complementary sensitivity function T(j@) can be achieved beyond this range. 

It is now assumed that also seeks a good performance (characterized by S(j@)) residing in the 
disturbance rejection. This performance can be achieved by finding the controller K(s) by solving the 
following equation: 


|| Ws(ieo).SGiw) |] 4 Vo €[0, ] (22) 

This equation 1s equivalent to: 

| WS(j@).SQG@) | œ < yconst Y wE[0, æf (23) 
Yconst 1S an arbitrary constant which is not necessarily equal to y (22). 
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Performance objectives and the robust stability can be simultaneously achieved by seeking a 
controller K(s) satisfying the following inequality: 


Maaa T 


W-0o).T0o) 








The choice of the weighting functions is as follows: Depending on the required performance, 
Ws(j@) is large at low frequency and smaller at high frequency, representing a low pass filter. 

In our case, it is a diagonal matrix Ws(j@) = ws.I, where wS represents the weighting function 
selected such that: 


1 
S(jo) ( iw, Go) | 
S or Oo (S(j — 
o (S(ja@)) < W, Go) 


for multivariable systems. 

Generally, the uncertainties and dynamics are neglected high frequencies, Wr0œ@) must be 
represented by a high pass filter to ensure robustness for high frequencies (œ > œ@0) and acceptable 
performance for low frequencies (@ < œ0). 

Our choice of weighting coefficients is completely connected to the parameter y as explained above 
(following the objectives set such that the bandwidth in this case equal 70rd/s [6]), after trial and error, 
a value of y = 1.92 was selected and the transfers functions of the weighting coefficients are: 


Ws = 10s+1 and W, 278 


~ 10(/70s+1) ~ 0.01s+0.01 


These weighting coefficients WS and WT define the controller K(s). The responses of the active and 
reactive powers are illustrated in Figure 10. From these results, it can be noted that the powers track their 
respective references perfectly, reflecting the right choice of parameters of our K(s). In addition, through the 
weights Wr and Ws defined, the interaction between the powers is completely eliminated and hence the 
controller was able to decouple the system and ensure good performance and achieve the control objective 
[18-26]. 
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Figure 10. Response of the active power (a) and reactive power (b) with K(s) control 


3.3. Robustness Test 
3.3.1. Increase of the Inductance of the Line by 20%: 

We have introduced changes to one of the line parameters (inductance) to analyze the robustness of 
our controller. Note that the increase of 20% of this parameter, does not affect the characteristic as shown in 
Figure 11 (a), because the active power follows perfectly the set point which explains the qualities of our 
controller K (s) Hoo. 
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Figure 11. Active power (a) and reactive power (b) corrected by the regulator Ho for the inductor of the 
line increased by 20% 


Figure 11 (b) shows the curve of the corrected reactive power after reduction of the inductance of 
the line by 20%, it is found that our characteristic is similar to that of Figure 10 (b), that is to say that it has 
not undergone any variation, which explains the robustness of our corrector K(s). 


3.3.2. Reduction of the Inductance of the Line by 20%: 

Figures 12 (a) and (b) represent respectively the curves of the active and reactive powers corrected 
after the reduction of the inductance of the line of 20%, these characteristics prove once more that our 
regulator retains its qualities, in particular its robustness. 
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Figure 12. Active power (a) and reactive power (b) corrected by the regulator Hoo for the inductance of the 
reduced line by 20% 


We can then conclude that the various changes made did not affect the qualities of the controller, 
which implies that the controller K(s) defined by the Hoo approach is perfect for our system. 


4. CONCLUSION 

In this work, we focused to show the control of the UPFC by the Hoo control. The control Hoo could 
ensure the desired performance, this implies that the control Hoo is robust and also implies a good synthesis 
of coefficients and weighting functions thereof. We can therefore probably say that the Hoo approach is, 
in the field of robust control, the theme that continues to give rise to a large number of publications following 
very great efforts. All this 1s not due to chance, because it represents a very natural way of formulating the 
problem of robustness. 
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